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1 Introduction 



A recent fit [Ij of Einstein-Cartan's tlieory [21 |3l H] to tlie Hubble diagram of supernovae 
was ratlier encouraging in tliat tlie spin-density, wliicli is tlie source of torsion, could 
- within today's error-bars - replace cold dark matter. The underlying space was fiat, 
homogeneous, isotropic and invariant under inversion. In the present work we extend our 
fit by admitting parity violation. Our motivation is that weak forces do break parity. 
Parity odd Einstein-Cartan theory in the context of maximally symmetric cosmology has 
a long history. Already in 1978, Bloomer |5] analyzed the theory, with the parity odd part 
only, on the 3-sphere. In 1986, Peter Minkowski [B] reconsidered it on fiat M"^. He mentions 
the handedness of spiral galaxies as motivation. More recently, in 2002 Capozziello et al. 
[7] took up the fiat, parity odd theory again. They claim that spin-density can replace 
the cosmological constant. 



2 Notations and field equations 

We use the conventions of reference [1] . For the reader's convenience we briefly summarize 
them. 

Let {fi = 0, 1, 2, 3) be a coordinate system on an open subset of M"^. We will use 
both a holonomic frame dx^ and an orthonormal frame, un repere mobile using Cartan's 
words, =: e°'^dx^, a = 0,1,2,3. We denote a metric connection with respect to an 
orthonormal frame by uj°'b ='■ oj°'bc^'^- It is as a 1-form with values in the Lie algebra of 
the Lorentz group. We follow the traditional convention and denote the same connection 
with respect to the holonomic frame by a different letter: F"^ =: F^^^dx'^, a gC,{^) 
valued 1-form. The link between the components of the connection with respect to the 
holonomic frame F and with respect to the orthonormal frame u is given by the GL{4) 
gauge transformation with e{x) = e^^^x) G (jL(4); 

oj = eFe"^ + ede~\ (1) 

Then (suppressing all wedge symbols) Cartan's two structure equations read: 

R:=duj + ^[uj,uj], (2) 

for the curvature 2-form i?"^ =: ^R'^hcd^'^^'^, and 

T := De = de + ue, (3) 

for the torsion 2-form, T° =: ^T°'bcG^e'^- It will be useful to decompose the torsion tensor 
into its three irreducible parts: 

Tabc = Aabc + VabVc - r]acVb + Mabc, (4) 

with the completely antisymmetric part Aabc '■= |(^afec + Tcab + Tbca), the vector part 
Vc '■= ^TabcV"'^y and the mixed part Mabc characterized by Mabc = —Macb, MabcV""^ = 0, 

and Mabc + M^ab + Mbca = 0. 
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In a Riemann-Cartan space there are two kinds of geodesies: curves that minimize the 
arc-length (or proper time) with respect to the metric and curves whose tangent vectors 
are parallel with respect to the connection. These two geodesies coincide if and only if 
the torsion only has a completely antisymmetric part, V — 0, M — 0. 

With these notations, the Einstein-Hilbert action reads 

where 60123 = 1- The energy-momentum current is the vector- valued 3-form Tq obtained 
by varying the orthonormal frame in the matter Lagrangian: 

Luie + f,uj]- jCuhuj] =: -/V^ + 0(f ). (6) 

The energy- momentum tensor Tab is defined by *Ta ='■ Tab^'^- 

Likewise, the spin current is the Lorentz-valued 3-form Sab obtained by varying the 
connection in the matter Lagrangian: 

CM[e,uj + x] - Cuie,^] =: -\r'Sa, + O(x'). (7) 

The spin tensor Sabc is defined by *Sab ='■ Sabc^'^, where * is the Hodge star of the 
metric. Einstein's equations are obtained by varying the total action with respect to the 
orthonormal frame: 

{R"'' + |Ae"e'') e'^eabcd = -IQttGtc or equivalently Gab - ^Vab = STiGTab, (8) 

with the Einstein tensor Gab '■= R'^ach — ^R'^'^cdVab- 

Likewise Cartan's equations are derived by varying the total action with respect to 
the connection: 

T'^e'^^abcd = —SnGSab: (9) 

or equivalently: 

Acab + ^Var]bc - ^VbTjac + M^ab = SnGSabc- (10) 



3 Homogeneous and isotropic spaces 

The invariance of the metric tensor g,j,i,{x) — e"^(x) e^v{x) rjab under an infinitesimal dif- 
feomorphism ^ is expressed by the Killing equation: 

Likewise the vector field ^ preserves the connection if 

^ dx- dx~^ dxi- '^''^ dx- "'^^ dx^^dx- ~ ^ ' 



The most general Riemann-Cartan space invariant under 5*0(3) xM^ has the Robertson- 
Walker metric: e° = dt, = adx, = ady, = adz, with the scale factor a(t), a pos- 
itive function of cosmic time t. The non-vanishing components w^bc of the most general 
5*0(3) K M.^ invariant connection are: 

h J' 

U^ij = U^Qj = -Sij, U^jk = — €ijk, (13) 
Qj Ct 

with two additional functions b{t) and f{t). The first is parity even like the scale factor, 
the second is parity odd. 

The Riemann tensor has the following non-vanishing components: 

-R^iOA: = -R*00fc = — Sik, R'^ijk = "2 Cijk, (14) 

. f . \? - P 

R^jOk = — ^ijk, R^jke = 2 i^ik^ji " ^ii^jk)- (15) 

Ct ct 

The Einstein tensor has: 

G™ = 3^, G„ = -(2- + -^j*., (16) 

The torsion tensor has: 

T^oj = — - — Sij, Tijk ~ 2 — eijk- (17) 

The antisymmetric part has only space components, Aij^ = 2f /aeijk, the vector part has 
only a time component, Vq = {b — a') /a, and the mixed part vanishes, M = 0. This 
result agrees with the curvature and torsion found in references [H [9] for spacetimes with 
maximally symmetric 3-spaces. 



4 Equations of state and Friedmann equations 

The most general 50(3) x M'^-invariant energy-momentum tensor contains two function of 
time, the energy density p{t) and the pressure p{t) and one usually assumes an equation 
of state p{t) =: w p{t). 

Likewise the most general 50(3) x M^-invariant spin density has two functions of time 
s{t) and s{t) in the two irreducible components: 5*0^^ =: —s{t) 6jk and Sijk ='■ —s(t) e-ijk 
and we assume two equations of state: 

s{t)=:wsp{t), ~s{t) =:wsp{t). (18) 

Then the generalised Friedmann equations, i.e. the tt and the xx components of Einstein's 



4 



equations, and Cartan's equations read: 



52 _ f2 

3 ^ = A + SnGp, (19) 

b' 52 _ f2 

2 - + = A - 8nGp, (20) 

a 

3 = StcGwsP 21 

a 

f 

2L = STcGwsp. (22) 
a 

These equations agree with results in references [HI E] and in reference [5] except for a 
missing 1/3 in front of the last term on the right-hand side of equation (23) there. Note 
that this factor re-appears correctly in the subsequent equation (24). 

However we disagree with a result in reference [6] stating that the field equations imply 
that the function f{t) (in our notations) must be constant. This result is reproduced in 
reference [7], probably because of a missing factor 3 in its equation (20) (in the arXiv 
version) and it is only with this factor missing that the presumably constant / can be 
interpreted as a cosmological constant. 

Putting the pressure to zero, p = 0, we have four equations for four unknown functions: 



a, b, f and p. Equations (19) and (22) are algebraic, the other two equations, (20) and 
(21), are first order differential equations for a and b. We use the algebraic ones to 
eliminate p and /. Then we have a unique solution with two inital conditions a(0) = Oq 
and fe(O) = bo. We therefore have five parameters, oq, 60, A, Ws and Wg. (We assume 



Newton's constant known.) These five parameters then fix p{0) by use of equation (19): 



2 9 2 

1 = ^mO + ^AO + 2f2sO ~ ^sQ + ^ ^sO> (23) 

with familiar dimensionless quantities: 

^ 8nGp ^ A ^ rr^^Gp ^ rr^^Gp , , 

■■= , := 7^ , ■■= WsH , VL-s := WsH . (24) 

In particular, we see that the scale factor today ao has dropped out. This is well-known for 
cosmology with vanishing spatial curvature and remains true in presence of non-vanishing 
torsion. Note also that the sign of / does not matter because only its square appears 
in the two Einstein equations. Therefore we may assume the state parameter ws to be 
non-negative. 



5 Hubble diagram 

To compute the Hubble diagram, we must solve the geodesic equations for co-moving 
galaxies and for photons [lO] . For both, torsion decouples and they reduce to the geodesic 
equations with the Christoffel connection of the metric. Consequently the redshift is still 
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given hy z = ao/ait) — 1 and the apparent luminosity d is still related to the absolute 
luminosity of the standard candle L by 



We have put the earth at the origin of the Cartesian coordinates and the supernova on 
the X-axis: 

xit):= / (26) 
Jt a{t) 

Note that the Einstein equations in presence of half-integer spin do feel torsion. How- 
ever the link between the Hubble constant Hq and d{z'^£)/dz{0) is purely kinematical and 
therefore does not depend on torsion. This fact will be crucial to identify consistently the 
initial conditions of Friedmann's equations. 



6 Data analysis 

The data analysis used in this paper has been fully described in [T] . Only a brief reminder 
is given here. The type la supernovae Hubble diagram is constructed using the Union 2 
sample [11] with 557 supernovae and the full systematic error matrix. The magnitude of 
supernovae is written as M{z) = ms + 2.51og£(z) where rris is a normalization parameter 



fitted to the data and £(z) the apparent luminosity (25). 



The apparent luminosity is computed using the generalized Friedmann equations (19) 



(20), (21) and (22). These equations are solved numerically by using the Runge-Kutta 
algorithm [T2] . 

The MINUIT [H] package is used to fit the best cosmology by minimizing the 
defined as: 

= AM^y-^AM, (27) 

where AM is the vector of differences between measured and expected magnitudes and 
V the full covariance matrix including systematic errors. 

Marginalization over unwanted parameters as rUg and error estimates or contour con- 
structions are obtained using the frequentist prescription [15]. The Einstein-Cartan cos- 
mology fit is performed with 3 or 4 free parameters (m^, Qm, ^s) while Q\ is derived 



from the Friedmann-like equation (23). 



Table 1 presents the results of the fit of Einstein-Cartan's theory (parity even and/or 
odd) and, for comparison, the results of the fit of the pure Einstein theory in the first line. 
Because of very high non-Gaussianity, errors are given at 1 and 2 sigma level. Minimum 

for all theories are statistically equivalent. If in the parity even case the preferred value 
for is compatible at a level of one sigma with baryon matter density, in the odd-parity 
case, the preferred value of Qm is in agreement with the total matter density of 0.27 
published by the WMAP collaboration [13]. This is not surprising since the preferred 
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Einstein 


— 7i — 1 l\ 1/1 1 (1 ic — 

n QK+U.iU+U.ib 
'-'•'-'^-0.11-0.17 


n QQ-t-u.iy+U.zo 
'-'•""-0.11-0.32 


0. 


0. 


530.0 


even-parity torsion 


'J-'J^-0. 07-0.08 


n OO+0.10+0.12 
'-'•"'-'-0.16-0.23 


n n/i +0.01+0.02 

'-'•'-'^-0.07-0.12 


0. 


530.4 


odd-parity torsion 


n 97+0.03+0.06 
'-'•^'-0.02-0.27 


n 70+0.04+0.06 
'-'• ' 'J-0.11-0.32 


0. 


n +0.22+0.55 
'-'•-0.22-0.55 


530.4 


odd-even parity 


n 08+0.^7+0.9 
'-'•'-'"-0.07-0.08 


n op:+0.10+0.15 
'-'•"^-0.15-0.25 


n 0/1+0.02+0.06 
'-'•'-'^-0.06-0.34 


n +0.01+0.6 
'-'•-0.01-0.6 


530.0 


odd parity no A 


n ni +0.02+0.03 

'-'•'-'-'--0.02-0.03 


0. 


0. 


n fifi+o.oi+0.02 

'-'•'J '-'-0.01-0.02 


560.7 



Table 1: Fit results (1 and 2a errors) for Einstein and Einstein-Cartan theories with even 
and odd parity. No flatness constraint is imposed in the pure Einstein's theory. 



value of ilg is exactly zero implying a flat Universe in the pure Einstein theory. In all 
cases, the cosmological constant energy density is only slightly changed. 

In figure [T]^a) the result of the Hubble diagram fit with odd Einstein-Cartan theory 
is shown (upper curve) with data points and error bars. As in the case of parity even 
Einstein-Cartan theory, the agreement between fitted curve and data points is excellent. 
The lower curve shows the fit resulting from putting the cosmological constant to zero for 
parity odd torsion. The agreement with data points seems good and suggests that the 
cosmological constant can be replaced by parity odd torsion. 

To test this hypothesis quantitatively, we use the log likelihood ratio technique. The 
log likelihood ratio is defined as: 

R = 2Ln ^^^^^^^"^^'^^ ^ (28) 
Here Sup denotes the supremum of the likelihood function defined in term of the '■ 

(27r)("/2)|V|i/2'' ' ^^^^ 

n is the number of data points and V the full error matrix. Thus the log likelihood ratio 
reads simply: 

(30) 

The probability distribution of this test variable is approximately a distribution with 
degree of freedom equal to the difference between the degrees of freedom of both models, 
one in this case. 

The minimum for the null cosmological constant hypothesis is equal to 560.7 while 
the minimum for odd parity torsion with cosmological constant is 530.4 (Table 1). The 
p-value is found to be equal to 610"*^ corresponding to 5.4 a significance. 

Because systematic errors of supernovae intrinsic magnitude variations at high redshift 
(above 1) can be important, we check the null cosmological constant hypothesis using only 
supernovae at redshift below 1. The log likelihood ratio is found to be 25 leading to a 
p-value of 5.310"'' or a significance of 5.01 sigma. Thus, the null cosmological constant 
hypothesis within the parity odd Einstein-Cartan theory is ruled out at more than 5 a. 
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For completeness, we perform the same analysis using simultaneously even and odd 
parity torsion. The resulting fit is slightly better because of one more fitted parameter. 
The of the fit is equal to 560.1 leading to a p-value of 5.9610"^ corresponding to a 5.4(T 
significance. 



7 Conclusion 

We find that a fit of Einstein-Cartan's theory to the Hubble diagram is incompatible at 
5 a level with the replacement of the cosmological constant by torsion, parity preserving 
or not. We think that the contrary claim in reference [7J relies on a wrong coefficient in 
the field equations [6l [7] . 
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Figure 1: (a) Fit results using the Union 2 Supernovae sample. The red (upper) curve corre- 
sponds to the Einstein-Cartan 3-fit (m<j, $7^,) ^^s) while the green (lower) curve represents the 
2-fit assuming a vanishing cosmological constant, (b) 39%, 68% and 95% confidence level contour 
in the (Omi^s) plane for the Einstein-Cartan 3-fit. 
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